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1 Introduction 



In a pioneer work |T[], H. Bethe solved the XXX Heisenberg model. The Hamiltonian of the 
model is written as 

N N 

H (J, h) = [SjS]^^ + S]S]^^ + +2hJ2 S- + constatn, (1.1) 

i=i i=i 

where J is the coupling constant and h is the external field. The energy of the system is 
given by 

^ 2J 

Em = 2W + ^^^, (1.2) 

j=i j 

where M is the number of up-spins and the variables xj are required to satisfy the Bethe 
Ansatz (BA) equation, 

' xj + ^ Y = Yl ~ ^ (13) 
Xj — i J Xj — Xji — 2i 

It was shown |]1| that every eigenstate of this model corresponds to a solution of this simul- 
taneous equation. 

The XXX Heisenberg model is the quantum integrable system. Thermal Bethe Ansatz 
method (TBA) and Quantum Transfer Matrix method (QTM) are the well- 

known methods to derive thermodynamic quantities for quantum integrable systems. TBA 
was originated by Yang and Yang , who derived the thermodynamics of repulsive 5-function 
bose gas system. Then, on the assumption of the so-called string hypothesis, M.Takahashi |P 
applied TBA to the Heisenberg model. Recently, some analytical results of thermodynamic 
quantities using QTM were presented by A. Kliimper 0,|l3]- Many pieces of concept with 
respect to QTM had been introduced before then. The equivalence of the results using TBA 
and QTM was proved by several researchers [llll,p^,|9|,p!3|. 



Our purpose is to derive the free energy of this system independent of TBA and QTM. 
In this paper, we use a method which we call a direct method. We start from the string 
hypothesis. There is no other mathematical assumptions except the string hypothesis. In 
other words, using the direct method, a problem based on the string hypothesis is solved. On 
the other hand, TBA is a method to solve the problem using a definition of non-equilibrium 
state entropy. As a main result of this paper, the direct method justifies TBA. 



The direct method was used for the ^-function bose system |jT^,|T^,|T6|,0. We apply the 
method to the Heisenberg model here. We define z and /5 as 2; = e"^^'' and /3 = The 
direct method enables us to express the cluster expansion coefficients hn in arbitrary order 
explicitly, 



log {z^'I^Ty e 



PH{J,h) 



:i.4) 



n=l 



and, the free energy. 



i- log [Tr e-^^(^''^)l 



1.5) 



by use of a function which is a solution of a nonlinear integral equation. Thus, we can show 



that the expression of the free energy perfectly agrees with the results of TBA [[TT . 

The outline of this paper is the following. In §2, we derive an explicit expression of cluster 
expansion coefficients in all orders. In §3, we prove that the free energy can be expressed 
in terms of a solution of a nonlinear integral equation. The last section is devoted to the 
concluding remarks. Technical details of calculations are summarized in Appendices A-H. 



2 The cluster expansion 



The string hypothesis is formulated as follows. In the thermodynamic limit, the energy of 
the Heisenberg model (|1 . 1|) can be written as 

2Jn 



E{J,h) = 2hM+J2. „.2^ 

n.a \ ci/ 



Here, the variables are determined by 



x'^ — m 



N 



{m,P)^{n,a) 



n _ m\ 



where 



E rj.m. ( ) 



(x — (n + m) i) (x — 


n — m 


^) 


(x + (n + m) i) (x + 


n — m 





min(n.m) 

n 



X + {n + m — 2l)i 
X — {n + m — 21) i 



(2.1) 



(2.2) 



(2.3) 



Equation ( |2.2| ) is referred to as string ceter equation. Every eigenstate corresponds to a 
solution of these simultaneous equations. Taking a logarithm of both sides of (|2.2|), we 



obtain 



2Ntan-'(^] = 27rI2 + A„,„ (x^ - x^) , (2.4) 

where 

/ N / . min(n,m) , - 

A„,^(a;) = -2tanM ) -2(l + 5„„)tani( )+V4tani( ).(2.^ 

\n + m J Wn — m\ J \n + m — 21 J 



From now on, we call ( p^.4D modified Bethe Ansatz (mBA) equation. 

The string hypothesis indicates that the partition function is a sum of the Boltzmann 
weights with respect to integers(or half integers) which define the energy of the system, 

^N/2rj,^ ^-PH{J,h) ^ ^g-/3£;(J,h)-2/3/iM_ ^^2.6) 

When we take the thermodynamic limit, we can replace the summations over discrete vari- 
ables by integral. 

There are some points which we have to consider carefully. The Fermi statistics must be 
taken into account in this integral. In this case, the exclusion property leads to the conditions 
that x'^ 7^ and that a state specified by {xj^} and a state specified by {x'^} which satisfies 
x''^ = a^p(Qf) same, where p is a permutation. 

To express a concrete form of ( ^.7| ), we prepare 3 symbols. 

definition 1. 

Na denotes the number of elements in a set A. 
definition 2. 

O (A) denotes all the patterns of division of a set A. A pattern of division 
is represented as a set having elements each of which is a cluster. The cluster is 
one of pieces into which a set A is divided, and the cluster is regarded as a set. 
For example, for A = {a, b, c}, Q{A) means 

0({a,6,c}) = {{a},{6},{c}},{{a},{6,c}}, 

{{b}, {c, a}}, {{c}, {a, b}}, {{a, b, c}}. (2.8) 



On condition that the domain of G(n) is the natural number, n G N, G(n) is 
interpreted as 6({1, 2, ■ ■ ■ , n}). 

definition 3. 

(B) denotes all the patterns of division of a set B, when each element of 
-B is a set and each pattern of division satisfies the following condition; all the 
sets which are included in each cluster have the same number of elements. For 
example, we can show for o"i = {a, b}, a2 = {c}, = {d}, 

e{Wi,a,,a,}) = {K},{a2},K}},{{ai},{a2,a3}}. (2.9) 

Using the above symbols, the integral (|2.7| ) can be written explicitly in the following form. 



n=l 




dl 


/ 


dx 



'■ 6»ee(n 



n( 

.6»'GC 



(2.10) 



We explain functions, ||^|^^ and E {() used in this equation, where ( G Q{6) and 9 G 9(n). 
The Jacobian ||^|^^ is constructed as follows. First, we consider 6 uniquely determined 
by the condition ( G G(^). Second, we make correspondence between each element a of 6' 
and two suffixes (n, a) for the variables and introduced in the string hypothesis. In 
this case, the condition A^^. = n is required. Third, we consider the mBA equation for the 
variables {1^}, {^a} which are constructed out of 9. Fourth, we introduce xgi, Igi, where 
9' is an element of (. Here, we have constraining conditions; xg' = x^ and Igi = I^, where 
(n, a) corresponds to cr G 9'. Finally, ||f is the Jacobian between {xg'} and {Ig'} which 
satisfy respectively the mBA equation made from 9 and the constraining conditions made 
from (. Explicit expression of such Jacobian matrix is given in eq. ( |A.1|) . We show the above 
procedure by an example. For the case. 



C = {0i,92}, 9, 



{{1,2}}, 9, 



{{3}, {4}}, 



(2.11) 



we have 



9 = {{1,2}, {3}, {4}}, ( e Q{9) 
5 



2Artan-i (^^) = 27rJe, + 2A2,i {xe, - 



(2vr)^ 



9/ 



AN 2N 
N,c ^01 + 4 a;^^ + 1 

m 



+ 



2A^ / 2 6 



X 



-23^^(7 ^^ + 7 ^^)- (2-12) 

^2 



+ 1 V(xe2 - +1 {xe2-xei) +9 



The energy E {() is expressed in terms of {x^} made from 6, 



= j:t^^.- (2-13) 



For the case (|2.11| ), E {() is given by 



^(0 = ^^ + 2^^. (2.14) 

Note that the constraining conditions made from ( have been used. 

Now, we show how ( 2.101) is derived. We notice two important facts. One is that, we 



must exclude unphysical states from the sum, when the sum is represented as the integral in 
( p.7|) . The unphysical state means the state corresponding to {x^} where some two elements 
of {x^} coincides, x^ = x^. Those states should be excluded because the norm of the state 
constructed out of {x^} satisfying such a condition is zero. In eq.( |2.1U| ), 

Hi-lf^'-'iNe'-iy. (2.15) 

6»'GC 

are factors which exclude the unphysical states. In the expression (|2.1CI|) , we add unphysical 
states in some terms of series and subtract them in some other terms. Then, the sum of pos- 
itive terms and negative terms is zero with respect to unphysical states. This mathematical 
manipulation enables us to write the free energy explicitly like ( p.lO|) . Let us explain this 
manipulation for the case of ( in ( |2.11| ) in detail. The states corresponding to ( in ( p.ll| ) 



are unphysical states. We define (' = {{1, 2}}, {{3}}, {{4}}. (' corresponds to the results 
of a string center equation ( |2.2| ) which defines one two-string 0:1^2 and two one-strings X3 X4. 
On the other hand, ( = {{1, 2}}, {{3}, {4}} in ( p.ll| ) corresponds to the results of the string 



center equation which satisfy X3 = X4. Therefore, the results corresponding to ( are part 

6 



of the results corresponding to Here, in the sum with respect to elements in ( 2.1CI| ) 



both ( and (' are summed up. And, the sum of coordinates ( |2.15|) of terms corresponding to 
both ( and (' is zero. Therefore, unphysical states corresponding to ( do not affect the sum 
in ( p.lOp . The other important fact is that, we must divide the sum by the number of the 
symmetry. Here, the symmetry means that when {x^} is a solution of the mBA equation, 
{x'^} which satisfies the relation x'^ = x^^^^^ is also a solution of the mBA equation. In this 
case, it is necessary to divide the sum by the number of the symmetry because {x^} and 
{x'^} correspond to the same state. In eq.( p.lO|) , 



/n\ (2.16) 

is a factor for the correction of the sum which takes such symmetry into account. 

We have met similar correction factors in the analysis of the 5-function bose gas system. 
The reason for the appearance of the correction factors is common. A detailed proof for the 
^-function bose gas system case is given in §3.1 ]TB| . 

In order to express the Jacobian ||f |jy^ explicitly, we further define three symbols. 



definition 4. 



Mb denotes a number of elements in one of the sets. Here, all the elements 
of a set B are sets, and all sets in B have the same number of elements. Then, 
Mb is defined as the number of elements in one of the sets. 

Me = N,, (2.17) 

where aeO eC^Q {9') and 9' E Q (n). 

definition 5. 

A {A) denotes all the patterns of connection of a set A. Here, what we call a 
pattern of connection satisfies the following two conditions. 1) Any two elements 
of A are connected or not. In other words, there is no multiple connection. 2) 
There is no closed path in the connections. Then, a pattern of connection is 
represented as a set of elements each of which corresponds to a connection. Here, 
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an element corresponding to a connection is a set which consists of two elements 
connected by the connection. For example, we have 

A{{a,b,c}) = iD,{{b,c}},{{c,a}},{{a,b}}, 

{{a, b} , {a, c}} , {{b, a} , {b, c}} , {{c, a} , {c, b}} . (2.18) 

definition 6. 

G{[X,A]) is an element of 6(A), where A G A{A). In other words, ^([A,^]) 
is a pattern of division of A. The pattern of division 6 satisfies the following 
conditions; 1) any two connected elements are in the same set in 6', and 2) elements 
of 6 are larger in number than elements of any other 6' G 0(^4) satisfying 1). For 
example, we can show that 



G([{{a,6}},{a,6,c}]) = {{a,6},{c}}. 
By use of the above three symbols, the Jacobian can be expressed as 



(2.19) 





dl 

dx 






-1 




(2vr)^c 








Y\ -NeNeiKMeMe' i^e - xe') 








AeA(C) 


{e,e'}ex 



X n 

C'6G([A,C]) 



E 



2NgMeN 
Ml + xl 



(2.20) 



where Kn^m {x) = d^n,m{x) / dx . A proof of eq.( p.20|) is given in Appendix A . 

At last, we show the explicit expression of the cluster expansion. Taking a logarithm of 
both sides of ( |2.10|) and rewriting the r.h.s., we get 



log r^Tr e-^'-)] = E 

n=i eee(n) 





dl 


/ 


dx 



-mo 



Cee(0) 



N.,-1 



1)\ 



c7V,C 



Y\ dxe'. 



(2.21) 



A proof is given in [Appendix B| . Here, 1 17 l^^v ^^^^ order term of 1 1^ |^ ^ regarded as a 



dx I cN,C, 

polynomial with respect to the number of spins, N . From the expression ( |2.21| ), the cluster 
expansion coefficient 6„ is given by 



bn 



- y 



eG0(n) 



E 



Cee(e) L^'gC 



JJ(_l)iv.'-i(Ar,,_i)! 





dl 


/ 


dx 



e Y\ dx0 



(2.22) 



where the exphcit form of | ^| ^ is 



dx I c7V,C 





dl 

dx 






-1 




(27r)^^ 




\{Ne 


E 


J]^ -NeNe'KMg.Mf), {xe - xq 




cN,C 




AeA,(C) 


{e,0'}ex 



2NgMgN 



Ml + x2 



(2.23) 



The symbol Ac (A) is defined as folfows. 
definition 7. 

Ac (A) is a subset of A (A), where any element A G Ac {A) satisfies the condi- 
tion Ng[\,a] = 1- For example, we can show that 

Ac({a,6,c}) = {{a,b},{a,c}},{{b,a},{b,c}},{{c,a},{c,b}}. (2.24) 



A derivation of the expression ||;;:|^^^ in ( p.23| ) is given in [Appendix C 

3 The free energy 

From eqs.( |2.1(]| ) and ( p.21| ), we can easily show that 



^1/2 [Tr 



'l3H{J,h)-\ 1/^ 



n=l 




dl 


/ 


dx 



E 



Ui-lf^'-HNo'-l] 

9'e( 



^ mo 



(3.1) 



where ||f a means |ff in which N is replaced with 1. 



,C I dx\N,C 

Now, we introduce a function u ix) 



u{x) = 1 + ^ z'^Un (x) 



(3.2) 



n=l 



where Un{x) is given by 



Un {x) 



- y 



6»ee(n) 



E 



J](-l)^«'-i(iV,,-l)! 

.6"eC 





dl 


/ 


dx 



\x] e 



-mo 



W dxe' 



6»'ec 



(3.3) 
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with the modified Jacobian | |f | ^ ^ (y) defined as 
dl 



dx 



(y) 



n 

CeG([A,c]) 



AeA(C) 



E 



{e,e'}e\ 
2NoMf> 



^ Ml + {xe - yy 



(3.4) 



.dec 

Comparing this definition with (|2.20|) , it is readily seen that ||f 1^ ^ = | |f ^ (0). Then, it is 
clear that we can get a relation 



l3H(J,h)-\ 



u{0)z 



-1/2 



(3.5) 



It can be shown (see [Appendix D| for a proof) that the function m„(x) satisfies the 
following recursion relation, 

2J/3 \ ^^^^ ( _ 2Jf3 



Un+1 (x) 



0+ 



expi-(^^j ^exp^-(— 



^ n! 



X — y — 2i ' X — y + 2i 



.(Tee 



dy 



Ui (x) 



exp 



2J/3 



exp 



X — y — 2i 



+ 



2J/3 



X — y + 2i 



dy 

2m 



(3.6) 
(3.7) 



where n is the natural number. Using this recursion relation we obtain a relation 



u (x) 



-2+1 + 



0+ 



exp 



2J/3 



(y-i)2+l 



X — y — 2i 



X — y + 2i 



z dy 
u iy) 2TTi 



(3. 



These results, ( |3.5| ) and (p.8[), are the same as those of |TT[] in the XXX Heisenberg model 
limit. 



4 Conclusion 

We have shown a method, which we call direct method, to derive the free energy of the XXX 
Heisenberg model using the BA equation. The cluster expansion coefficient in arbitrary order 
is given for the first time. The expression of the free energy perfectly agrees with that by 
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TBA. It is remarkable that the free energy is obtained without recourse to TBA and QTM. 
On the other hand, there remains a problem. We have started from the string hypothesis 
in this paper. Therefore, it is still a challenging problem to derive each cluster expansion 
coefficient only from the BA equation. This problem is under investigation. 

We think, however, that these results have theoretical significances. What we have done 
in this paper is that we calculate a certain summation using combinatorial argument with 
mathematical justification. The summation may contain a problem which comes from the 
string hypothesis. On the contrary, TBA solve this problem as follows. We define a non- 
equilibrium entropy, and minimize the free energy using such entropy. Then, the condition 
of minimization gives the summation. In other words, TBA entirely relies on the physical 
definition of the entropy. In this paper, it is proved that the two results by TBA and by 
the direct method are the same. Therefore, in addition to the (5-function bose gas, for the 
XXX Heisenberg model, the direct method gives the mathematical justification of TBA. To 
summarize, it is suggested that TBA can be systematically proved for many other integrable 
systems using combinatorial argument presented here. 
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Appendix A Explicit form of Jacobian 



dl 



dx 



In this appendix, we prove ( p.20|) . Recall that ||f is the Jacobian defined for ( G 6(^), 
where 6 G Q{n). The Jacobian matrix can be rewritten as follows; the Jacobian matrix has 
the order x N(^, and each element of the matrix is 



27r 



dxf) 



2MeN AT 



M. 



[xe - xe" 



{x) = dAn,m{x)/dx. 



where 9, 6', 6" G C and Kr, 

Here, we generalize the setting of the problem. We define N x N matrix 

en - Zlfc^n Ckan,k if n = 771 
Cnan,m if U ^ 171, 



{A.l] 



(A.2) 



where {ctn^m}, {c„} and {e„} are arbitrary series and a„_m is symmetric, a„ m = c^m.n- This 
matrix is a generalized form of the Jacobian matrix (|A.1|). Therefore, what we have to prove 



becomes that the determinant of A„ ^ is 



■ N 

n 

n=l 



n -1 



E 

AeA({l,-,Af}) 



{n,m}£A 



n 

aeG([A,{l,-,7V}]) 



E 



.(A.3) 



First, we prove that each term of the determinant |v4„^m| can be considered as a pattern of 
connection A G A ({1, ■ ■ ■ , N}). When we regard the determinant as a multi- variable polyno- 
mial with respect to {a„_m} and regard connection between n and m, we can consider 
each term of the polynomial as a pattern of connection A. In other words, we can prove the 
following two facts. First, there is no multiple connection. That is to say, there is no an,m to 
the power of 2 or more in any term of the polynomial. Second, there is no closed path made 
from connections. That is to say, there is no product ap(i),p(2)a.p(2),p(3) ■ ■ ■ '3'p(M-i),p(A/)ap(M),p(i) 
in any term, where p is a permutation. 

To prove these two facts, we use the method of false position. We assume there is a 
term that has at least one closed path through M connections. Without loss of generality, 
the above assumption can be rewritten that there is a product 01,202,3 ■ ■ ■ om-i.a/Oa/,! in any 
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term. Then, we introduce a matrix A'^ ^ such that 

- Efe^n ^kan,k \in = m<M 
en - Zlfe^n Ckan,k H u ^ m > M 



A' = 

n.m 







a n<M<m or m<M<n 
other. 



(A.4) 



Now, we regard ^ multi-variable polynomial with respect to {an,m}- From the defi- 

nition of determinant, it is easily shown that the term with a product 01,2^2,3 ■ ■ ■ cim-i,m^m,i 
in has the same coefficient as a term with such product in However, 

is identically because the matrix A'^,^ is linearly dependent. Therefore, there is no term 
containing a product 01,202,3 • • • om-i,mOm,i- Substituting 2 for M, we can prove that there 
is no Qn^m to the power of 2 or more in any term of the polynomial. Then, we can consider 
each term of the matrix as a pattern of connection A. 

Second, we find where a term considered as a pattern of connection is formed. From the 
definition, we can write the determinant as 



N 



sgnp • Yl 

pePN 



(A.5) 



n=l 



where Pjv is a permutation group with respect to {1, • • • , N}. It can be shown that any term 
considered as a pattern of connection is in an expansion of a product only in the case p is the 
identity permutation. The following is a proof of the fact by the method of false position. If 

there is such term in the case p is not the identity permutation, there exist no and Iq which 
satisfy p'°(no) = no, p^(no) ^ uq and ^ 0, where /q > / and p^in) means that p operates 
on n for m times, pip{- ■ - pin) ■■■))■ Therefore, we can rewrite 11^=1 ^n,p{n) as 



7o-l 

S'("o)^pH"o),p'+H"o) 

.1=0 



(A.6) 



This term clearly has a closed path, and confiict with the definition of the p. 
Finally, we can rewrite |v4„ „i| as follows. We expand 

N 
n=l 



(A.7) 
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and regard it as a multi- variable polynomial with respect to {an,m}- |^n,m| is a sum of terms 
which are in this polynomial and can be considered as a pattern of connection. This means 
that the term has no multiple connection and no closed path. This sum means the r.h.s. of 

Appendix B A proof of ( |2.21| ) 



In this appendix, we prove ( p.21| ) from ( |2.1CI| ). 
We define by 



n=l 



From (|2l0| ) we get 



-y 

ege(n) 



E 

Cee{e) 



JJ(_l)A^«'-i(iV,,-i)! 

6»'ec 









dl 


I 


dx 



(B.i: 



where ||f is the Jacobian defined for (. It is convenient to introduce f(Q as 



/(C) = 





dl 


/ 


dx 



mo 



(B.3) 



Using this function and the relation (|C.1|) proved in [Appendix C| , we can rewrite Zn as 

= ^ E 

n! ^-^ 

6»G0(n) 

We also define 6„ by 



Tee 



Cee(e) 



JJ(_l)A^«'-i(iV^,_l)! 

6»'eC 



E n/(o 



(B.4) 



log 



n=l 



E"' 



n=l 



2; . 



(B.5) 



We can easily show the following relation between 6„ and Zn, 

= A E iNo-m-lf^-'l[NJZ 



eee(n) 

Substituting ( |B.4D for Z„ in ( [B.6| ), we get 



(B.6) 



6»Ge{n) 
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n 



9'g0(<T) 



Cee(e') 



JJ(_l)A^«"-i(iV,„-l)! 



- y 



9Ge(n) 



.cre6» 



E 

Cee(e) 



- y 



HA'. 



E 



.e'GC 



E 

«6e(c) 



E n/K') 

5ee{C) C'e? 



n / (CO 



/(C). 



(B.7) 



In the first equality, we have only performed the substitution. In the second equality, we 
have done the following manipulation. We regard each side of the equality as the sum with 
respect to elements in a set which satisfies some conditions. Therefore, we may change the 
order of the sums. Here, we show a simple example of "change of order". 
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5^5^/({n,m}) = ^5^/({n,m}). (B.8) 

n=l m=l m=l n=m 

Similar to this example, both sides of the second equality are the sum of the same function 
with respect to elements in the same set. In the third equality, we have simplified the last 
term of the l.h.s. using a relation. 
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Finally, combining (|B.1|) , (p.3| ), ( p.5| ) and (p.7| ) we arrive at the equation (|2.21 
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Appendix C The structure of the Jacobian 

In this appendix, we prove 
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and eq. (|2.23|) from eq.( |2.20|) . We here repeat that ||f |^^^ is the first order term of ||^|^^ 
with respect to A^. 

Eq. (|2.2CI| ) can be rewritten as 
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In this equality, we have changed the order of the sums. We can easily show the relation 
( p.23| ) from this expression. Then, eq. (|C.2|) and eq. (|2.23|) prove 



Appendix D A proof of (|S^) and (jS^) 

In this appendix, using the mathematical induction we prove the relations with respect to 
M„(x), ( ^.61 ) and (|3.7]). The definition of 'U„(x) is given in eq.( 
We can change the r.h.s of ( |3.7D into 

2 f 2JI3 \ dy 



l + ix-yf"""^^ V y' + iy 2vri" 
This is the same as the definition of Ui{x). Thus, we have shown the relation ( p.7| ). 

We suppose ( |3.6| ) and ( |3.7| ) hold for 1 < n < riQ, where hq is a natural number. For 
simplicity, tlq is merely written as n in what follows. Then, the function Um {x) has poles 
only at X = ±2i, and satisfies lim|^|^oo = 5m,i, where 1 < m < n. By use of these two 

properties, the r.h.s. of (|3.6|) becomes 
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In this equality we have changed the path of integration. Both C+ and C_ are the paths of 
integration which counterclockwise surround the region '^y > 0. When these paths are in 
the neighborhood of real axis, C+ passes above the real axis and C„ passes below the real 
axis. We separate ( |D.2| ) into two parts, and calculate them separately. 
First, we consider the first term of ( |D.2| ) which is rewritten as 

Ne\ r ^ sY{^^,-N,\uNAy-si) -^dy 
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The integral diverges in case that we do not sum up with respect to s. However, the sum 
turns out not to play an important role in the following calculation. Hence, to shorten the 
expression, a symbol X]s=±i omitted for a while. Now, by use of the relation ( [E.l|) proved 
in [Appendix E| , eq. ( p.3|) is rewritten as 
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On condition that we take a limit 5^0 and Mg = 1,1/ (M| + (a;^ — y + si)'^) yields a pole 
on the real axis with respect to the variable y. To detour the pole on the real axis and 
symmetrize the integral paths, we use the relation (|F.1|) proved in Appendix F . Then, ( D.4 ) 
becomes 
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Here, we assume the relation + 1} G cr„+i G 9n+i G C- We also use this symbol in the same 
sense as in ( p.l2| ), ( p.l9| ) and ( p.20| ). The constant i?„ in (p.5| ) is known as the Bernoulli 
number and is determined by 
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For example, first few -B„'s are 
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The symbol J^^^ in ( p.5|) indicates that each path of integration symmetrically avoids poles 
on the real axis. Because all the poles on the real axis have the form 1/ {xg — xqi), we can 
write explicitly f as 
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The functions C {6, 6') and E {6, 6') are defined by 
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Note here that the second term of (p.lO|) is for E {6, 9') used in ( p.5|) because of the relation 
^6»„+i = 1- Then, by use of the relation 
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which is easily shown, ( p.5|) becomes 



- y 



9ge(n) 



Tee 



E 

Ce0(e+{n+i}) 



TT(_i)«.-.W--l)'1 ^ / 

A^fl' Ne 

g/^^ o' J J sym 



n 

.6»'eC 



JJ C {6', 6" 
9',e"}ex,e',e"^0n+i 



{e„+i,e'}e\ 



. (D.12) 



Changing the order of the sums, this expression becomes 
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Thus, using the above relations we have shown that (p.3|) is expressed as ( P.13 ). 
Second, we write the sum of the second and third terms of (|D.2D , 
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where, C±i means C±. In the following change of the expressions, the sum with respect to 
s does not play an important role. This situation is similar to the one for ( p.3|) . Then, 
the symbol Xls=±i ^^^^ omitted for a while. Using the relation ( PD.II) , we have from 
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We use a relation 
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, 27r ^_ + ^3^)2 

where /o-(x) and /(x) are analytic functions in the regions —1 < Qx < and < Qx < m 
in case of s = 1 (0 < Qx < 1 and —m < Qx < in case of s = —1) respectively, and 
maXo-ge ^o- = fn- Then, ( p.l5| ) becomes 
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We point out that in this case the first term of the r.h.s. in ( p. 16 ) is because there is no 
pole inside (outside) the path of integration C+ (C_) with respect to y including the point 
\y\ = oo. Then, using a simple relation J2eee„ Ylaee (^)^'^~^ (^o- — 1)' = 5n,i, we get from 
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K' + (a^eo - + s {Me, -l)if _ 
where = Wo}- Here, we have changed the path of integration from Cg to —oo + s5i 
oo + s6i, because there is no pole at \xg,\ oo. In eg. ( |D.18| ), when we take a limit 5^0 
and Mg = Mg, + 1, the last term contains a pole on the real axis with respect to x^q. 
Symmetrizing by use of ( |F.1|) the integral paths which avoid the pole on the real axis, we 
obtain from ( p.l8|) . 
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Changing the order of the summations, this expression becomes 
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From the above relations we see that ( p.l4| ) is equal to ( p.21| ). 

Both parts of the r.h.s. of ( p.2|) are replaced with (p.l2|) and (p.21|) . Then, the r.h.s. of 
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Using eq. (|G.l|) proved in [Appendix Q , we can confirm the definition of Un+i{x)^ which is the 
r.h.s. of dU). 



Appendix E An expression of Un{x) 

In this appendix, we prove the following expression of Un{x), 
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Here, the definition of u„ (x) is ( p.3| ), and the definition of ||f [^^^ ^ (y) is 
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Substituting ( p.3| ) for UNa{y) in the l.h.s. of ( p. 1|) and using a relation 
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which is proved in the same way as (C]T), we can show that the l.h.s. of ( |E.1| ) is calculated 

as 
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In this equality, we have changed the order of the sums. Using a relation 
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for simplification of the first term of ( IE .41 ), we obtain the r.h.s. of ( |E.1| ). 

Appendix F Symmetrization of the paths of integrals 



In this appendix, we prove a relation. 



dxfo (x) Yl 



n=l 



dyn 1 
2m yn-x 



fn iVn) 



E 

(7C{1,---,N} 



f /o (x) 


n (^) 




n 


J sym 






ne{l,---,N} 



1 . ( ^ 

In yVn) „ . 

yn- X 2m 



dx. (F.l) 



Here, the definitions of i?„ and J^^^ are in ( |D.6| ) and (|D.8| ) respectively, and we assume that 



fn{x) is an analytic function on the real axis. 

The paths of integrations in the r.h.s. of ( [F.l| ) are changed into ones in the l.h.s. of 
( [F.l| ). In other words, symmetrical paths of integrations are changed into ones which keep 
the relation < Then, the r.h.s. of (|F.1|) becomes 
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Using the following two relations which are easily shown, 
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for (|F3), we obtain the the l.h.s. of (FU|) . 
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Appendix G Alternative representation of the modi- 
fied Jacobian 



In this appendix, we show that 
dl 



dx 



(y) 



.dec 



y y y - 

^ ^ ^ sMe - [xe - v) i 

AeAc(C) eec s=i-i y \ y y) 



9',6»"}ga e'ft"^e 



n ^(^'^') 

{e,e'}e\ 



{G.l] 



is an alternative representation of the modified Jacobian ( |3.4| ). Here, C {0,6') and E {6,6') 
are defined in (p.9|) and ( P.10|) respectively. 
We can rewrite ( |3.4| ) as 
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In the first equality, we have changed the order of the sums. In the second equality, we have 
used the relation ( |H. 1|) proved in [Appendix H| . Applying a relation 
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to ( |G.2| ), we obtain the r.h.s. of (|G.1|) . In the relation ( |G.3|) , it is supposed that s is ±1. 
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Appendix H A Proof of a fractional relation 



In this appendix, we prove 
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where a„, c„, x„ and y are arbitrary numbers 
We first prove a relation 
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We regard the r.h.s. as a polynomial of degree — 1 with respect to y. It is clear that this 
relation holds on the A^ points y = Xn + SnO.n'i- Therefore, this equation is nothing but an 
identical equation. Dividing both sides of this identical equation by Hnt'^n"'^™ (-^n — u) i]dn^, 
we get 
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Here, we restrict that s„ is ±1. The sum of each side of this equation with respect to {s„,} 
gives 
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We are in a position to prove the main issue of this Appendix. Using ( |H.4|) for the Lh.s. 



( [H.l| ), we get 
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Expanding the r.h.s. of ( |H.1|) , we also get (pi .51 ). This completes a proof of (|H.1| ). 
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